NATIONAL  ADVISORY  COMMITTEE 
FOR  AERONAUTICS 


TECHNICAL  NOTE  2378 


AUTOMATIC  CONTROL  SYSTEMS  SATISFYING  CERTAIN 

GENERAL  CRITERIONS  ON  TRANSIENT  BEHAVIOR 

By  Aaron  S.  Boksenbom  and  Richard  Hood 

Lewis  Flight  Propulsion  Laboratory 
Cleveland,  Ohio 


Washington 

June  1951 


DUO  qualhy  4 


Reproduced  From 
Best  Available  Copy 


20000816  163 


UISTRIBlJtfON  STATEWE^?T  A 
Approved  for  Public  Release 
Distribution  Unlimited  jt 


NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 


TECHNICAL  NOTE  2378 


AUTOMATIC  CONTROL  SYSTEMS  SATISFYING  CERTAIN 
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By  Aaron  S.  BoksenBom  and  Richard  Hood 


SUMMARY 

An  analytic  method  for  the  design  of  automatic  controls  has 
“been  developed  that  starts  from  certain  arbitrary  criterions  on 
the  behavior  of  the  controlled  system  and  gives  those  physically 
realizable  equations  that  the  control  system  can  follow  in  order 
to  realize  this  behavior •  The  criterions  used  are  developed  in 
the  form  of  certain  time  integrals. 

General  results  are  shorn  for  systems  of  second  order  and  of 
any  number  of  degrees  of  freedom.  Detailed  examples  for  several 
cases  in  the  control  of  a  turbojet  engine  are  presented. 


INTRODUCTION 

In  the  past  several  years,  there  has  been  increasing  develop¬ 
ment  and  interest  in  automatic  control^  in  the  fields  of  gun 
direction,  guided-missile  control,  and  control  of  gas-turbine 
engines,  for  instance,  where  very  refined  and  accurate  controlled 
behavior  is  required,  need  still  exists  for  further  development  of 
the  methods  of  controls  analysis  and  design. 

Recent  developments  in  this  field  have  been  mainly  concerned 
with  the  problem*  of  control  analysis  both  in  the  realm  of  linear 
systems  (reference  l)  and  in  the  realm  of  nonlinear  systems  (ref¬ 
erence  2).  These  analytical  works  answer  the  following  questions: 
How  will  a  given  system  behave  or  how  is  its  behavior  affected 
when  certain  constants  in  the  system  are  changed? 

,  Another  problem  of  equal  and,  in  some  cases,  greater  import¬ 
ance  is  that  of  control  synthesis.  Work  on  this  problem  seeks  the 
answer  to  the  following  question:  Given  certain  criterions  con¬ 
cerning  the  behavior  of  a  controlled  system,  how  should  the  system 
be  designed?  The  answer  to  this  question  should  give  all  aspects 
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of  "the  system}  for  instance^  whether  the  system  should  he  linear, 
what  the  general  configuration  should  he,  and  what  the  precise 
values  of  all  the  constants  should  he. 

This  sjmthesis  problem  has  hardly  heen  broached  in  literature. 
The  use  of  analysis  as  a  design  procedure  offers  a  partial  solu¬ 
tion  to  this  problem  in  that  the  analysis  of  a  large  number  of  cases 
may  reveal,  coincidently,  a  case  that  satisfies  the  desired  cri- 
terions  of  controlled  behavior.  Such  a  method  is,  at  best,  long 
and  tedious  and  almost  always  would  result  in  compromises  because 
the  systems  chosen  to  be  analyzed  wotild  probably  be  such  that  they 
could  never  satisfy  all  the  desired  criterions. 

A  method  for  designing  a  linear  system  to  satisfy  certain 
special  criterions  when  operating  on  a  random  input  is  developed 
in  reference  3.  This  method  is  applicable  as  an  addition  to  a  con¬ 
trol  system  whenever  random  external  disturbances  are  involved. 

A  partial  solution  to  the  synthesis  problem  is  developed  in  refer¬ 
ence  4  in  satisfying  the  criterions  of  noninteraction  for  systems 
with  many  degrees  of  Treedora. 

An  analysis  made  at  the  NACA  Lewis  laboratory  and  presented 
herein  develops  a  rational  method  of  control  synthesis  that  starts 
from  any  arbitrary  but  physically  realizable  criterions  and  results 
in  the  equations  for  the  best  system  that  satisfies  these  criter¬ 
ions,  As  is  shown,  the  nature  of  the  criterions,  in  general, 
requires  minimizing  certain  time  integrals  by  using  the  calculus 
of  Variations  and  the  methods  developed  are  an  application  of  the 
calculus  of  variations  to  the  problem  of  control  synthesis. 

A  careful  scrutiny  is  first  made  of  the  whole  problem,  fol¬ 
lowed  by  a  development  of  general  results.  These  general  results 
are  then  applied  in  examples  to  the  design  of  turbo jet -engine  con¬ 
trol  systems.  In  general,  the  methods  used  vary  according  to  the 
order  of  the  differential  equation  describing  the  plant,  the  proc¬ 
ess,  or  the  system  being  controlled  and  according  to  the  nxmiber  of 
degrees  of  freedom  or  independent  vairiables .  Detailed  analyses 
are  presented  for  application  to  a  system  of  first  order  and  of 
one  degree  of  freedom.  General  results  for  systems  of  second 
order  and  any  manber  of  degrees  of  freedom  are  developed  in  the 
appendix. 
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SYMBOLS 


The  following  synibols  are  used  in  this  report: 


a,b,c 

C 

E 

F,H 

Fy(y,y) 

F(ty) 

fj^^fg,  .  .  . 


gas -turbine -engine  characteristic  constants 

constant 

function  of  'K 

functions  of  y  and  y 

indicates  partial  differentiation  with  respect  to  y 
indicates  F  [y(ti),  ^(ti)] 
arbitrary  function 

weighting  functions  used  for  gas-turbine  control 
criterions 


G,G^,G2 


functions  of  y 


g 


function  used  for  gas -turbine -engine  surge  criterion 


h 


function  used  for  gas-turbine-engine  blow-out  criterion 


K 


controller  gain 


L,Li,L2 


temperature  limits  for  gas-turbine  engine 


N,P,T,w^ 


N 

e 


deviations  of  gas -turbine -engine  speed,  compressor- 
discharge  pressure,  characteristic  temperature, 
and  fuel  flow,  respectively,  from  values  at  some 
common  equilibrium  condition 

actual  engine  acceleration 


n 


exponent 


t 


time 


bti 

v,y:,z. 


time  at  end  of  transient 

variation  in  time  at  end  of  transient 

independent  variables,  functions  of  time 
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€ 

a 


T 

Subscripts : 


small  number 
arbitrary  constants 

gas-turbine-engine  time  constant  for  response  to 
teit5>erature 

transient  time  constant  of  controlled  system 


initial  condition  of  variable 

final  value  of  variable 

setting  or  desired  valtie  of  variable 


Superscripts ; 

*  indicates  case  different  from  optimum 

The  dot  indicates  differentiation  with  respect  to  time. 

The  prime  indicates  differentiation  with  respect  to  the  argu¬ 
ment  shown. 


ANALYSIS 

Survey  of  Problem 


Control  problem.  -  An  important  aspect  of  the  control-synthesis 
problem  is  a  clear  definition  of  the  criterions  of  desired  con¬ 
trolled  behavior.  If  a  variable  y  is  to  be  controlled,  a  reason¬ 
able  criterion  is  that  the  time  integreil  of  some  function  of  y 
is  to  be  a  minimum  or  a  constant;  that  is. 


pti 

I  f(y)  dt  =  constant  or  minimum  (l) 

UO 
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or 


nti 


(y-yg)^ 


dt  =  constant  or  ininimum 


(2) 


Equation  (2),  for  instance,  weights  the  error  in  y  as  the  square 
and  according  to  the  time  duration  of  that  error.  Another  type  of 
criterion  may  he  that  which  requires  a  certain  time  duration  to  he 
a  minimum  or  a  constant;  that  is, 


Ati 


u 


dt  =  constant  or  minimum 


(3) 


The  use  of  a  single  criterion,  such  as  equation  (l),  will  usu¬ 
ally  yield  f(y)  =  constant.  This  result  is  reasonable  because 
f(y)  can  usually  he  made  identically  a  constant  if  no  additional 
criterions  are  imposed  on  other  variables  in  the  system.  Usually, 
certain  limiting  conditions  exist,  however,  on  other  variables  in 
the  system  and  these  conditions  must  be  included  in  the  original 
criterions . 


Thus,  a  possible  criterion  could  be  written  as  follows: 


for 


•Yq)  dt  =  minimum 


f(z)  dt  =  constant 


}  (4) 


If,  for  instance,  y  =  engine  speed  and  z  =  characteristic  tem¬ 
perature  of  a  gas-turbine  engine,  the  criterion  of  equation  (4) 
states  that  it  is  desired  to  design  a  control  system  such  that, 
for  a  particular  value  of  a  temperature  integral,  the  integral  of 
the  speed-error  squared  is  a  minimum.  This  criterion  may  be  used 
if,  for  instance,  it  is  knovn  that  an  over-ten^Jerature  condition 
can  be  tolerated  for  a  certain  period  of  time  and  it  is  desired  to 
keep  the  average  speed  error  at  a  minimum  during  a  transient. 
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The  general  theory  vill  shov  that  as  many  criterions  as 
desired  of  the  type  shown  in  equations  (l)  to  (4)  can  "be  included 
together  and  a  control  system  can  "be  derived  that  automatically 
satisfies  all  these  criterions  simultaneously. 

Another  aspect  of  the  control  criterions  is  the  end  conditions 
of  the  integrals  of  equations  (l)  to  (4).  The  time  interval  for 
which  these  integrals  are  to  be  a  minimum  or  a  constant  must  be 
chosen.  A  reasonable  time  interval  is  any  duration  during  which 
essential  external  disturbances  are  constant  and  during  which  the 
system  to  be  controlled  moves  from  one  essential  level  of  opera¬ 
tion  to  another.  The  essential  external  disturbances  are  those 
that  cannot  be  inmediately  corrected  by  the  control  system.  If  an 
essential  external  disturbance  were  allowed  within  the  time  inter¬ 
val  of  the  criterions,  no  physically  realizable  system  could  be 
expected  to  anticipate  this  disturbance  so  as  to  behave  properly 
before  this  disturbance  occurs.  An  essential  level  of  operation 
is  any  specific  condition  of  only  those  variables  that  must  be 
continuous.  It  will  be  shown  that  the  essential  level  of  opera¬ 
tion  appears  an  a  natural  boundary  condition  for  the  type  of  cri¬ 
terion  used.  In  the  case  of  a  turbojet  engine  the  tr^sient 
behavior  of  which  can  be  described  by  a  first-order  differential 
eauation,  the  engine  speed  determines  the  level  of  operation.  If 
a  "lag  exists  in  the  fuel  system  or  between  temperature  ^d  engine 
speed,  then  both  engine  speed  and  acceleration  are  required  to  des¬ 
cribe  the  essential  operating  level  of  the  engine . 

Analytic  problem.  -  The  control  system  resulting  from  any 
design  method  kist  be  physically  realizable.  There  are  two  aspects 
to  this  problem.  First,  it  is  possible  to  set  down  criterions  that 
are  not  realizable  with  any  system  or  are  incompatible  with  each 
other.  If  such  criterions  are  used,  the  -unrealizability  will 
appear  either  as  a  requirement  on  the  control  to  look  ahead  into 
the  future  or  bis  an  inability  to  satisfy  the  boundary  conditions 
of  some  differential  equation.  In  most  cases,  a  clear  understand¬ 
ing  of  the  criterions  used  and  of  the  system  to  be  controlled  will 
indicate  incompatibilities  of  this  sort. 

The  second  aspect  of  physical  realizability  is  purely  mathe¬ 
matical.  It  is  desired  to  derive  a  description  (a  differential 
equation)  of  the  control  or  the  controlled  system  that  satisfies 
the  criterions  of  control  and  all  the  necessary  boundary  conditions 
that  arise  in  the  derivation  of  this  equation.  Although  the  mathe¬ 
matical  solution  of  the  problem  may  be  any  derivative  or  integral 
of  this  differential  equation,  the  physical  solution  of  the  prob¬ 
lem  requires  the  differential  equation  that  itself  satisfies  the 
boundary  conditions  and  for  which  no  Tindetermined  constants  of 
integration  exist.  Thus,  such  forms  as 
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are  not  necessarily  interchangeable  as  descriptions  of  some  part  of 
a  controlled  system  because  the  forms  differ  by  an  undetermined 
constant  of  integration.  For  stable  linear  systems,  the  effect  of 
this  constant  becomes  vanishingly  small;  for  the  general  nonlinear 
systems  presented  herein,  however,  this  constant  must  be  considered. 

Stability  problem.  -  The  requirement  of  stability  is  a  special 
criterion  that  does  not  enter  into  the  main  body  of  the  methods  of 
this  report.  It  may  enter  in  the  final  steps  of  the  method  where 
the  final  differential  equation  describing  the  controlled  system 
may  be  the  integreil  of  a  higher -order  differential  equation  that 
satisfies  the  necessary  boundary  condition  for  stability.  In  addi¬ 
tion,  it  is  always  necessary  to  add  to  the  controlled  system  a 
stability  device  that  does  not  affect  the  behavior  of  the  system 
as  far  as  satisfying  the  other  cri ter ions  is  concerned.  This 
device  can  be  described  as  follows: 


when 

y  = 

1 

then 

y  =  0 

>  (6) 

or,  for  a 

second-order  system. 

when 

y  =  y3 

N 

and 

li 

o 

/ 

)  (7) 

then 

o 

II 

General  Theory  and  Results 

It  has  been  shown  that  the  criterions  for  control  can  be 
developed  in  the  following  forms: 
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rvti 


u 


f(y)  dt 


rvt- 


(y-yfi)  <i'fc 


)  (8) 


and  so  forth.  If,  for  such  a  list  of  criterions,  one  of  the  inte¬ 
grals  is  to  be  a  minimum  under  the  condition  that  the  other  integrals 
are  to  be  constant,  it  is  sufficient  (reference  5)  to  make 


The  A's  are  arbitrary  constants  that  enter  into  the  control  sys¬ 
tem  as  the  adjustable  parameters  and  are  precisely  determined  by 
the  choice  of  values  that  the  constant  integrals  are  to  have. 


The  technique  of  the  Tv  multipliers  is  widely  used  for  prob¬ 
lems  of  this  type  where  one  condition  is  to  be  a  minimum  under 
other  restrictive  conditions.  Indeed,  the  conditions  need  not  be 
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in  integral  form  and  any  functional  or  differential  relation  among 
variables  can  be  handled  in  a  similar  manner  (reference  6) . 

Equation  (lO)  can  be  made  very  general  -when  all  possible 
restrictive  conditions  are  included.  In  the  final  equations,  which 
are  derived  later,  if  any  one  criterion  is  not  to  be  used,  then  the 
corresponding  A— ►O.  If  any  of  the  criterions  are  to  be  zero  (as 
for  the  case  of  a  variable  having  an  absolute  limit),  then  the  cor- 
re  spending  A  — > 


If,  for  this  development,  the  system  to  be  controlled  is  of 
first  order  and  of  one  degree  of  freedom  (has  one  independent  ^vari¬ 
able),  then  the  variables  y  and  z  are  such  that  z  =  z(y,y). 
Equation  (lO)  can  then  be  written,  in  general,  as 


1 


^*(7^7)  =  minimum 


(11) 


In  equation  (U),  F  is  a  continuous  function  of  y  and  y,  and 
y  is  a  continuous  function  of  time. 

The  calculus  of  variations  (reference  5)  is  used  to  determine 
y  as  a  function  of  time  such  that  the  integral  of#  equation  ( 11 ) 
is  a  minimum;  that  is,  if  the  solid  curve  of  figure  1  makes  the 
integral  a  minimum,  any  other  curve  (such  as  the  dashed  one)  will 
make  the  integral  equal  to  or  greater  than  the  first  integral.  If 
the  two  curves  are  very  close,  the  condition  that  the  two  inte¬ 
grals  are  equal  will  make  the  integral  of  equation  (ll)  stationary 
(maximum,  minimum,  or  inflection  point). 

The  problem  is  similar  to  that  of  finding  the  maximum  or  mini¬ 
mum  point  on  a  curve  by  setting  the  first  derivative  equal  to  zero. 
VJhether  a  maximum  or  minimum  point  exists  is  decided  by  the  second 
derivative  at  that  point.  In  the  variational  problem,  proving  a 
true  minimum  involves  taking  the  "second  variation."  As  in  the 
problem  of  finding  a  minimum  on  a  curve,  the  second  variation  proves, 
at  most,  a  local  minimum  and  not  an  absolute  minimiun.  For  the 
specific  examples  discussed  herein,  an  absolute  minimum  is  proved 
by  another  method.  In  many  cases,  the  physical  meaning  of  the 
equations  will  indicate  the  existence  of  a  unique  minimum  that 
obviates  going  further  than  the  "first  variation." 


V 
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On  figxire  1,  the  curve  for  minimizing  the  integral  of  equa¬ 
tion  (ll)  is  the  solid  curve  having  the  value  y(t)  at  any  time. 
Any  other  curve  differing  from  it  hy  a  small  amoimt  is  shown  as  the 
dashed  curve  having  the  value  y(t)  +  €  &y(t)  at  any  time,  where 
c  is  a  small  number  and  &y(t)  is  an  arbitrary  continuous  fxmc- 
tion  of  time.  The  condition  for  the  integral  of  equation  (ll)  to 
be  stationary  is 


5t]_ 


- 

‘^*Jt=0 


F(y+e  6y,  y+c  5y)  dt  =  0 


(12) 


at 


e=  0 


vhere  — [5y(t)3  =  5y(t)  (reference  5) .  The  time  duration  for 
dt 

the  integral  of  equation  (12)  is  such  as  to  start  at  some  definite 
time  (t=0)  but  not  to  end  at  a  definite  time  but  rather  along 
some  curve  y  =  f(t)  (fig.  l),  in  order  to  allow  the  proper  botuid- 
ary  conditions  of  moving  from  one  essential  level  of  operation  to 
another,  as  previouslj^  discussed. 

Performing  the  operation  indicated  by  equation  (12)  leads  to 


pti  pti 

J  Fy  &y  dt  +  I  Fy  6y  dt  +  F(ti)  St^  =  0  (13) 

Integrating  the  second  term  of  equation  (13)  by  parts  gives 

^1 

+  F(tjL)  btj^  =  0  (14) 

0 


Because  Sy  is  an  arbitrary  function,  the  Integral  and  the 
boundary- condition  terms  must  vanish  separately.  The  geometry  at 
the  end  condition  (t=ti,  fig.  l)  gives 
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&y(ti)  =  |f'(ti)  -  y(ti)J  8ti  (15) 

as  e— >0.  Ihe  two  conditions  that  follow  from  equation  (14)  then 
"become 


and 


pti 


(16) 


F-(ti: 


[f'(ti)  -  y(ti)j 


-  Fy(0)  5y(0)  =  0 

(17) 


The  time  interval  during  which  the  criterion  of  equation  (ll) 
is  to  hold  is  considered  as  that  during  which  the  system  moves  from 
one  essential  operating  level  to  another j  in  this  case,  from  one 
definite  vsdue  of  y  to  another  definite  value  of  y.  Thus, 


5y(0)  =  0 


and 

f’(ti)  =  0 


(18) 


Equation  (l6)  is  satisfied  only  if  the  integrand  is  zero^  and, 
hecause  btj^  /  0,  the  two  conditions  of  equations  (16)  and  (17) 
hecoine 


Fy  =  ^  (yp 

and 

F(ti)  =  Fy(t3^)  y(ti)  (191) 

Equation  (l9a)  need  not  hold  at  t  =  0  because  &y(o)  =  0,  The 
only  condition  that  need  hold  at  t  =  0  is  that  Fy(o)  is  finite, 

which  will  be  true  if  y  is  continuous  at  t  =  0.  At  the  start 
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of  a  new  transient,  y,  F,  F  ,  and  F-  may  be  discontinuous, 

whereas  at  other  points  (O  <  t  <  t^^)  Fy  will  be  continuous 
because  of  equation  (19a). 

Equation  (l9a)  is  the  differential  equation  for  the  y(t) 
that  satisfies  the  original  criterion  of  equation  (ll) .  The  physi¬ 
cal  answer  to  the  problem  is  the  first  integral  of  equation  (l9a), 
which  satisfies  the  boundary-condition  equation  (l9b) .  This  solu¬ 
tion  is 


F(y,y)  =  y  Fy(y^y)  (2o) 

and  whenever  y,  Fy,  and  so  forth,  are  continuous. 


by  differentiating  equation  (20)  \d.th  respect  to  time.  Thus, 
either  y  =  0  (which  is  only  true  during  stability)  or  equa¬ 
tion  (19a)  is  satisfied. 

Thus,  equation  (20)  is  the  description  of  that  physically 
realizable  system  the  behavior  of  which  will  automatically  and 
simultaneously  satisfy  those  criterions  included  in  the  function  F 
during  that  time  interval  for  which  the  external,  disturbances  are 
constant  and  during  which  the  system  goes  from  one  operating  level 
to  any  other  operating  level.  A  stability  device  must  be  added  to 
the  system)  the  description  of  such  an  ideal  device  is  as  follows; 

■vdien 


y  =  yg 

then 

y  =  0 

An  additional  condition  must  be  met  if  y  is  discontinuous  in  the 
interval  0  <t  <  tj^.  In  this  case,  F^  must  be  continuous  during 

the  discontinuity  in  y;  Fy,  however,  will  usually  be  discontinuous 
when  y  is  discontinuous.  This  discontinuity  in  Fy  usually 

means  that  some  essential  external  disturbance  has  entered  the  sys¬ 
tem  and  the  point  of  discontinuity  must  be  the  start  of  a  new  time 
interval . 


NACA  TN  2378 


13 


Application  to  Turbojet -Engine  Controls 


In  the  usual  case  of  designing  turbojet -engine  controls,  the 
engine  speed  that  sets  the  essential  operating  level  of  the  engine 
and,  in  the  main,  other  pertinent  characteristics  such  as  thrust, 
is  to  be  set  or  controlled.  Limiting  conditions  of  the  engine  are 
those  of  overspeed,  overt emperatxire,  compressor  surge,  and  rich 
burner  blow-out.  The  following  criterions  on  the  behavior  of  this 
engine  are  typical: 

^t-, 


\‘tl 


\ti 


f^(N-Ng)  dtj  for  speed  control 


f«(N)  dt;  for  speed  overshoot 


f^CT)  dt;  for  temperature  overshoot  and 
undershoot 


f4[P  -  g(N)]  dt;  for  compressor  surge 


f5[p  -  li(N)]  dt^  for  blow-out 


dt;  for  rise  time 


>  (23) 


Figures  2(a)  to  2(e)  show  the  natiire  of  the  functions  appeairing  in 
equation  (23).  The  variable  involved  is  essentially  weighted  by 
the  function  shown  and  integrated  with  respect  to  time.  The  quan¬ 
tity  P  -  g(K)  is  the  amount  the  compressor-discharge  pressure 
exceeds  the  safe  pressure  for  surge  and  g(N)  (shown  in  fig.  2(f)) 
is  the  coirpressor-discharge  pressure  for  each  engine  speed  at  a 
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safe  value  belov  svirge.  Rich  burner  blow-out  can  be  handled  in  a 
similar  manner  (h(tl)  is  shown  in  fig.  2(g)).  The  rise  time  is 
the  total  time  for  the  system  to  move  from  one  essential  operating 
level  to  the  other. 


The  linearized  engine  characteristics  can  be  expressed,  by 
assuming  first-order  behavior,  as  follows: 

T  =  aN  +  affN  (24a) 

P  =  bN  +  cT  (24b) 


Thus,  the  integrals  of  equation  (23)  are  all  of  the  form 

\ti 

f(N,N)  dt 


(25) 


where  f  is  a  continuous  function  of  N  and  N,  and  N  (barring 
impulsive  temperatures  and  the  like)  is  a  continuous  function  of 
time. 


Speed  control;  case  A.  -  If  only  the  error  in  speed  control  is 
considered  important,  the  criterion  becomes 


I  f2^(H-Ng)  dt  =  minimum 

Jo 


(26) 


where  F  =  f2(N-¥s).  Equation  (20)  becomes 

f^(R-Ng)  =  0 

and  from  the  nat\ire  of  f^^  (fig.  2(a)), 


(27) 


This  result  means  that,  in  the  absence  of  other  criterions  on  the 
engine  behavior,  this  speed  control  should  keep  speed  error  identi¬ 
cally  zero,  vhich  is  physically  realizable  only  in  the  sense  of 
allowing  infinite  temperatures  and  the  like  to  keep  the  speed  error 
identically  zero.  This  result,  however,  is  inconsistent  with  the 
previous  development  of  equation  (20)  in  that  N  is  now  a 
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discontinuous  fimction  of  time  and  the  time  interval  of  the  inte¬ 
gral  of  equation  (26)  is  zero.  This  instance  is  actually  a  trivial 
case  of  the  general  problem.  The  result  (equation  (27))  does  indi¬ 
cate  that  a  criterion  like  that  of  equation  (26)  must  be  accompanied 
by  an  additional  criterion  (temperature  overshoot,  for  instance) 
to  give  a  physically  realizable  system. 


Speed  control  with  temperature-limiting  criterions;  case  B,  - 
If  the  error  in  speed  control  and  the  overshoots  and  undershoots 
in  temperature  are  to  be  considered  as  the  primary  criterions  of 
control,  then  from  equations  (lO)  and  (ll),  r 


pti 


u 


0 


Af,(T) 


] 


dt  =  minimxm 


(28) 


where  F  =  f^(N-Ng)  +  A  f2(T)-  Equation  (20)  then  becomes 

fj^(N-Ng)  +  A  f3(T)  =  A  aaN  f3’(T)  (29) 

and  the  ideal  stability  device  is  such  that 
when 

N  =  Ng 

then 

N  =  0 

Equations  (29)  and  (30)  describe  the  complete  control  system. 

In  figure  2(c),  it  is  convenient  to  let  f3(T)  =  (T-Lg)^  for 

T  >  Lg  and  f2(T)  =  (T-L^)^  for  T  <  L^.  In  general,  the  power  n 

should  be  >1,  because,  when  n  !^1,  T  may  be  infinite  and  of  such 
nature  as  to  make  N  discontinuous  and  physically  unreal  even 

r 

though  ^  finite.  In  the  examples  of  this  report, 

n  =  2  and  f^(N-Ng)  =  (N-Ng)^. 


% 
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A  variety  of  methods  of  setting  up  the  control  system  to  real' 
ize  equation  (29)  exist.  When  the  preceding  expressions  for  f^ 

f3  are  used^  equation  (29)  can  he  put  in  the  convenient  form 


(N-Ng)  ,  ^2  2  212 

- - -  +  (L-aN)  =  a  n  N 


where ^  for  acceleration^  when  N  <  Wg^  then 

9 

^  N  ^  0  and  L  = 

and^  for  deceleration^  when  N  >  then 

9 

U  <  0  and  1  =  1]^ 


/  (31) 


A  schematic  diaigrara  of  the  system  is  shown  ^  in  figure  3,  where 
equation  (5l)  is  considered  to  give  a  desired  N  for  value 
of  N  ,  and  L.  Consistent  with  equation  (3l),  a  right -triangle 

construction  is  used^to  give  a  desired  N.  The  actual  N  can  he 
made  very  close  to  N  hy  using  a  high-gain  proportional  controller, 
as  shown.  Provision  must  he  made  to  change  the  si^  of  N  and 
the  value  of  L  when  N-Ng  changes  sign.  In  addition,  the  sta¬ 
bility  condition  requires  a  provision  for  making  ^  =  0  whenever 
speed  error  is  very  small  or  zero. 


T'he  use  of  a  high-gain. proportional  controller,  which  follows 
from  the  requirement  that  N  he  discontinuous,  me^s  that  the 

fuel-flow  rate  required  may  he  infinite  if  lags  exist  in  the  fuel 
system  or  in  the  feed-hacks.  But,  as  no  criterion  has  been  set  on 
fuel-flow  rate,  this  requirement  does  not  violate  the  original 
criterions.  If  necessary,  however,  a  criterion  on  fuel-flow  rate 
may  then  he  added  to  equation  (28) .  Even  thou^  a  criterion  on 
ten^erature  is  being  satisfied,  figure  3  does  not  use  any  direct 
measurement  of  teicperature .  Actually,  the  equation  for  tempera¬ 
ture  (equation  (24a))  is  used  as  an  indirect  indication  of 
temperature- 

The  control  system  of  figure  3  has  one  adjustable  parameter  \ . 
For  any  value  of  this  system  will,  for  the  value  of  integral 
temperature -overshoot  squared  obtained,  give  the  minimum  value  of 
integral  speed-error  squared.  The  value  of  Tv  determines  the 
actual  value  of  the  integral  temperature -overshoot  squared. 
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The  integral  temperature-overshoot  squared  as  a  function  of  Tv 
is  shown  in  figure  4  for  the  special  case  where  aWg  =  Lj  that  is, 
acceleration  or  deceleration  to  the  speed  that  corresponds  to  lim¬ 
iting  temperature.  In  this  case,  the  system  of  figure  3  becomes  a 
simple  first-order  lag  system  and  equation  (31)  becomes 

E(L-aN)  =  aaii  (32) 


where 


In  figure  4,  the  integral  speed-error  squared,  the  maxiimnn 
ten^erature,  and  the  time  constant  for  this  transient  are  also 
shown  as  functions  of  A.  A  curve  showing  the  minimum  rise  time 
for  the  corresponding  temperature  integral  is  included  for  compari¬ 
son  that  will  he  discussed  later.  The  equations  for  the  curves  of 
figure  4  are  as  follows:  ^ 


(time  constant)  ^ 


I  j 


The  left  side  of  these  equations  have  been  put  in  dimensionless 
form.  The  maximiun  temperature  occurs  at  the  beginning  of 

the  transient.  The  time  constant  is  that  for  the  controlled 
system  and  is  shown  compared  with  the  engine-time  constant  a. 

From  figure  4,  the  value  of  A  is  chosen  as  a  compromise 
between  the  various  quantities  of  equation  (33).  For  E  =  l(^=soo), 
the  temperature  does  not  overshoot,  the  speed  integral  is  0.5,  and 
T  =  a.  As  E  increases  (A  decreasing),  the  temperature  integral 
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and  the  maximum  temperature  increase,  vhereas  the  speed  integral 

anH  the  time  constant  decrease.  A  value  of  E  =  (a  X  =  1) 

appears  to  he  a  good  compromise  and  is  used  for  the  subsequent 
discussion. 

The  behavior  of  the  system  of  figure  3  can  be  best  seen  by 
drawing  the  trajectories  in  the  phase  plane  shown  in  figure  5, 
where  aoN  is  plotted  against  aN  for  lines  of  constant  aUg 
according  to  equation  (3l) .  Lines  of  constant  temperature  are 
45°  parallel  lines  in  this  plot  and  the  lines  of  T  =  L  are  shown. 
Each  trajectory  intersects,  and  is  tangent  to,  the  line  T  =  L  at 
N  =  Ng.  Figure  5  completely  describes  the  transient  behavior  of 
the  system.  For  any  starting  point  anywhere  on  the  plane  (for 
instance,  point  A),  the  system  will  automatically  move  the  oper¬ 
ating  point  to  that  trajectory  corresponding  to  the  Ng  that 
exists  (point  B)  and  then  along  this  trajectory  to  the  point  T  =  L, 
(N  =  Hg,  point  C),  and  finally  the  stability  condition  will' enter 
to  move  the  operating  point  along  the  solid  vertical  line  to 
N  =  0  (point  D) . 

The  time  sense  for  these  transients  is  obtained  by  solving 
the  differential  equation  (equation  (31))  for  the  speed  and  the 
tenperature  time  responses .  The  equations  for  these  solutions  are 
as  follows : 


(34) 


These  transients  are  shown  in  the  general  case  in  figure  6  for  a 
step  increase  in  Ng .  Maximum  temperature  occurs  at  t  =  0  and 
the  temperatvire  overshoot  decreases  as  W  increases  and  when 
N  =  Ng,  T  =  L  and  T  =  0.  The  stability  condition  then  causes  T 
to  drop  to  its  equilibriiim  value  at  this  point.  The  time  scale 
shown  in  figure  6  corresponds  to  the  specific  relative  values  of 
the  ordinates  shown. 

Minimum  rise  time  and  temperattire-limiting  criterions;  case  C.  - 
In  order  to  obtain  a  minimum  rise  time  for  a  constant  temperature 
integral,  the  requirement  is  that 
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pti 

^  |l  +  'I'*'  =  minimum  (35) 

vhere  F  =  1  +  Xf3(T).  Equation  (20)  iDecomes 

1  +  X  fjCT)  =XaaN  fj'CT)  (36) 

and  the  ideal  stability  device  is  such  that 
when 

then 


where  Ng  is  the  desired  stable  N.  Equations  (36)  and  (37)  des¬ 
cribe  the  coEplete  control  system.  Equation  (36) ,  and  therefore 
the  control  system  proper^  does  not  include  Ng  and,  in  this  case, 
the  stability  device  is  independent  of  the  control  system  proper. 

V/hen  the  same  f3(T)  is  used  as  was  used  in  the  previous 
case,  equation  (36)  can  be  put  in  the  convenient  form 

1  /T  2„2-2 

+  (L-aN)  =  a  a  N 

where,  for  acceleration,  when  N  <  Ng,  then 

N  >  0  and  L  =  L2 

and^  for  deceleration,  when  N  >  Ng,  then 

* 

N  <  0  and  1  =  13^ 

A  schematic  diagram  of  this  system  satisfying  equation  (38)  is 
shown  in  figure  7.  This  control  system  is  the  same  as  the  previous 

system  (fig.  3)  except  that  in  figure  7  replaces 
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one  leg  of  the  ri^t-trlangle  construction  and  the  stability  con¬ 
dition  imist  be  imposed  outside  the  control  system,  as  Kg  does 
not  enter  into  the  criterion  of  control  (equation  (35)). 

The  system  of  figure  7  has  one  adjustable  parameter  %  which, 
as  before,  sets  the  precise  values  of  the  integrals  entering  into 
the  criterion  (equation  (35)),  as  well  as  all  other  behavior 
characteristics . 

The  temperatvire  integral  is  shown  in  figure  8  as  a  fuhction 
of  for  the  special  case  where  aWg  =  Lj  that  is,  for  accelera¬ 
tion  or  deceleration  to  the  speed  that  corresponds  to  limiting  tem¬ 
perature.  In  addition,  the  rise  time,  the  maximum  temperature, 
and  the  integral  speed-error  squared  are  shown  for  this  transient. 
A  curve  showing  the  minimum  speed  integral  for  the  corresponding 
temperatu.-e  integral  (from  case  B)  is  also  shown  and  vrill  be^dis- 
cussed  later.  The  equations  for  the  curves  of  figure  8  are  e.s 
follows : 


^-^0  /v/r(L-aKo) 


/\/L+7\(L-aKo)^ 

2(L-aKQ)VA’ 


- - - -  sinh”^  (L-aNQ) 

2X(L-aKQ)^ 


(39) 

Again  X  is  to  be  chosen  as  a  coi®)romise  between  these  various 

1 

quantities.  In  this  case,  has  the  tmits  of  temperature 

and  V^(L-allo)  is  the  pertinent  dimensionless  parameter.  For 
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a  fixed  7,^  the  initial  speed  Nq  vill  therefore  change  the  time 
integrals  and  an  increased  L-bUq  will  decrease  the  temperature 

integral  and  the  maximum  temperature  and  increase  the  rise  time 
and  the  speed  integral. 

The  dynamic  hehavior  of  the  system  of  figure  7  is  ^  shown  in  the 
phase  plane  in  figure  9.  This  figure  is  a  plot  of  auN  against  ail, 
according  to  equation  (58),  for  various  values  of  ^  .  Lines  of 
constant  temperature  are  45^  parallel  lines  and  the  lines  T  =  L 
are  shown.  For  any  fixed  only  two  trajectories  are  followed: 
one  for  accelerating  and  one  for  decelerating.  From  any  starting 
point  on  the  plane,  the  system  will  automatically  move  to  that 
trajectory  corresponding  to  acceleration  or  deceleration  'and  will 
move  along  this  trajectory  until  the  stability  condition  enters 
(at  N  =  Ng)  to  maJ^e  N  =  0. 

The  dependence  of  the  time  integrals  on  Nq  may  require  X 
to  vary  with  Nq,  In  figure  9,  Nq*s  corresponding  to  each  X  are 
shown  such  as  to  keep  X(L-aNo)^  =  8.  The  value  of  this  parameter 
was  taken  so  as  to  have  the  temperature  integral  in  this  case 
equal  to  that  of  the  previous  case  for  purposes  of  comparison. 


The  time  sense  for  these  transients  is  obtained  by  solving 
the  differential  equation  (38),  for  the  N-t  and  T-t  transients. 
The  equations  for  these  solutions  are  as  follows: 


Vl(ali-L)  = 
'\pX{T-L)  = 


(40) 


These  transients  are  shown  for  the  general  case  in  figure  10  for  a 
step  increase  in  Ng.  The  temperature  will  jump  to  some  value 
above  L,  b\it  in  this  case,  unlike  the  previous  case,  the  tempera¬ 
ture  continues  to  increase  as  the  speed  increases.  Maximum  tem¬ 
perature  occurs  at  the  end  of  the  transient,  ^^henever  N  =  Ng, 
two  such  conditions  are  shown,  the  stability  condition  takes  over 
and  T  is  reduced  to  its  equilibrium  value.  The  time  scale  shown 
on  figure  10  corresponds  to  the  specific  relative  values  of  the 
ordinates  as  shown. 

Comparison  of  cases  B  and  C.  -  In  case  B,  a  system  was  devised 
that,  for  a  constant  value  of  the  integral  temperature-overshoot 
squared,  gives  a  minimum  integral  speed-error  squared.  In  case  C,  a 
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system  >ra,s  devised  that,  for  a  constant  temperature  integral, 
gives  a  minimum  rise  time.  For  case  B,  the  rise  time  is  not  a  min¬ 
imum  hut  can  he  compared  vith  the  minimum  rise  txme  of  case  C . 

For  the  special  transient  of  accelerating  to  T  =  L,  the  system  of 
case  B  reduced  to  a  first-order  system  the  time  constant  of  which 
is  shown  in  figure  4.  Because  five  time  constants  are  considered 
as  the  rise  time  of  an  exponential,  the  time  constant  of  case  B 
should  he  compared  with  one-fifth  the  minimum  rise  time  of  case  C. 
In  figure  4,  for  corresponding  values  of  the  temperature  integral, 
one-fifth  the  minimum  rise  time  from  case  C  is  plotted.  Figure  4 
now  shows  that  the  rise  time  for  case  B  is  ahout  twice  the  minimum 
possible  rise  time. 


At  corresponding  values  of  the  temperature  integral. 


J(N-Nf )  dt  for  case  C  can  he  compared  with  the  minimum 
0 

possible  value  of  this  quantity  from  case  B.  The  minimum  value 
of  I  (N-N^)  dt  is  plotted  in  figure  8.  It  is  seen  that,  for 

J  0 

the  same  temperature  integral,  |  (N-N^)^  dt  when  minimum 


rise  "time  is  obtained  is  about  120  percent  of  the  minimum  possible 

pti 

value  of  I  (N-N^)^  dt. 


Proof  of  absolute  minimum.  -  As  previously  noted,  the  minimi- 
zation  of  equation  (11)  involves  not  only  the  condition  that  the 
first  variation  be  stationary,  which  leads  to  equation  (20),  but 
also  that  the  second  variation  be  positive.  This  second  condition, 
however,  would  prove,  at  most,  a  local  minimum  for  the  condition 
of  equation  (ll).  A  special  proof  of  an  absolute  minimum  is  shown 
as  follows : 

If  equation  (ll)  is  written  in  the  following  form: 

yf 

*  cLv 

F(y,y)  -r-  =  minimum 

y 


(41) 
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and  y  is  _ considered  as  a  function  of  y  for  the  minimtun  condi¬ 
tion  and  y*  is  considered  as  a  function  of  y  for  any  other  pos 
sihle  case,  then  the  condition  for  an  absolute  minimum  is  that 


F(y,H  ^  -  F(y,y)^>0 

y*  Uy^  y 


F(y,y*)  -  ^  F(y, 

L  y 


where  dy/y*  is  a  positive  differential.  If  the  function  F  is 
considered  in  a  general  form^  quadratic  in  y,  as  follows: 

F(y,y)  =  G(y)  +  y  Gi(y)  +  y^  GgCy)  (44 

then  equation  (43)  becomes 


(y*-y)  a2(y)  - 

uyq  L  y  J  y* 


The  equation  of  the  control  system  (equation  (20))  becomes, 
for  the  form  of  F  assumed  in  equation  (44), 

y^  G2(y)  =  G(y)  (46) 

Using  this  expression  for  G(y)  puts  equation  (45)  in  the  form 

pVf 

I  G2(y)  (y^^-y)^  %  >0  (47) 


For  cases  B  and  C,  G2  =  ?^a^cr^  where  0.  Thus  an  absolute  mini¬ 

mum  is  proved  for  these  two  cases. 
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Degree  of  minimum,  case  B.  -  The  left  side  of  eg.  nation  (47) 
can  he  used  as  a  measure  of  the  deviation  from  optimum  conditions 
when  the  methods  of  this  report  are  not  used.  If  the  two  cases 
(W  for  the  optimum  case  and  N*-  for  any  other  case)  are  compared 
at  the  same  value  of  the  temperature  integral,  the  left  side  of 
equation  (47)  is  the  difference  between  the  integral  speed-error 
squared  for  any  case  and  the  minimum  possible  value  of  this  quantity. 
The  ratio  of  this  deviation  to  the  minimum  value  becomes 


Fractional  increase  _  2E 
in  speed  integral  " 


IM 


<iy 


(48) 


for  the  'bransient  of  acceleration  or  deceleration  to  limiting  tem- 
perature.  The  coefficient  of  the  integral  of  equation  (48)  for 

the  value  of  X  previously  chosen  (a  X  =  l)  is  2^/2. 


Case  B  can  also  be  considered  as  giving  a  minimum  value  of 
the  temperature  integral  for  any  definite  value  of  the  speed  inte¬ 
gral.  If  two  cases  (N  for  the  optimum  and  N*  for  any  other 
case)  are  compared  at  the  same  value  of  the  speed  integral,  the 
left  side  of  equation  (45)  is  proportional  to  the  difference 
between  the  integral  temperature-overshoot  squared  for  any  case 
and  the  minimum  possible  value  of  this  quantity.  This  deviation 
can  be  written  as 


Fractional  increase  in  _ 
temperature  integral 


py'f 

2E  I  M-T  y  dy 
(E-l)^JyQ  Vf"'^C/ 


(49) 


for  the  transient  of  acceleration  or  deceleration  to  limiting  tem¬ 
perature.  The  coefficient  of  the  integral  of  equation  (49)  for 

the  value  of  X  previously  chosen  (a^  =  l)  is  16.4. 

Degree  of  minimum^  case  C .  -  If  the  two  cases  N  and  are 

compared  at  the  same  value  of  the  temperatxire  integral,  the  left 
side  of  equation  (47)  is  the  difference  between  the  rise  time  for 
any  case  and  its  minimum  possible  value.  This  deviation  can  be 
written  as 
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Fractional  increase 
in  rise  time 


■K  (L-aWo)^ 
s  ii]h~^  (L- sNq  ) 


ay 

ay* 


(50) 


for  the  transient  of  acceleration  or  deceleration  to  limiting  tem- 
peratrire.  The  coefficient  of  the  integral  of  equation  (50)  for 
the  value  of  /s/T(L-aKQ)  =  previously  chosen  is  4.54. 


Case  C  can  also  he  considered  as  giving  a  minimum  value  of 
the  temperature  ^integral  for  any  definite  value  of  rise  time.  If 
the  tvo  cases  N  and  N*  are  now  compared  at  the  same  value  of 
rise  time,  the  left  side  of  equation  (47)  is  proportional  to  the 
difference  between  the  integral  temperature-overshoot  squaired  for 
any  case  and  the  minimum  possible  value  of  this  quantity.  This 
deviation  can  be  written  as 


(51) 


for  the  transient  of  acceleration  or  deceleration  to  limiting  tem¬ 
perature.  The  coefficient  of  the  integral  of  equation  (51)  is  16.4 
for  the  value  of  'v/T(L-aNo)  previously  chosen. 


GENERAL  SUMMARY 

VJhen  the  criterions  on  the  behavior  of  a  controlled  system  can 
be  expressed  in  certain  general  forms,  as  follows: 


« 
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r»t 


0 


f(y)  dt 


pi 


)  (8) 


where  the  time  interval  is  taken  as  any  duration  during  which 
essential  external  disturbances  are  constant  and  during  which  the 
system  moves  from  one  essential  operating  point  to  another,  the 
optimum  control  can  be  considered  as  one  that  minimizes  one  of  the 
integraJ-s  of  equations  (8)  while  maintaining  the  other  integrals 
at  prescribed  values.  The  analytical  problem,  according  to  the 
calculus  of  variations,  reduces  to  the  following  equation: 


r*i 


f(y)  dt  +  (y-ys)^  dt  +?V2  I  ^0^2)  dt  +  A3  I  dt  =  minimum 


0 


0 


(9) 


For  general  first-order  systems,  where  z  =  z(y,y),  equation  (9) 
reduces  to 


F(y,y)  dt  =  minimTim 


(11) 


where  F  =  f(y)  +  IviCy-yg)  H-Ag  fo(z)  +  A3.  The  equation  neces. 
sary  for  the  control  system  to  satisfy  equation  (U)  and  all  the 
boundary  conditions  becomes 

F(y,y)  =  y  F‘(y,y) 


» 


(20) 
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This  equation  should  he  followed  hy  the  control  system  proper.  In 
addition,  a  stability  device  must  he  added  to  the  system  the  ideal¬ 
ized  characteristics  of  which  would  make  y  =  0  when  y  =  yg. 

The  arbitrary  multipliers  X  are  then  found  hy  evaluating 
the  integral  criterions  involved  in  F.  The  transient  behavior  of 
the  system  derived  is  found  hy  solving  the  differential  equa¬ 
tion  (20).  The  degree  of  the  minimum  or  the  amount  suffered  when 
any  other  control  system  is  used  was  evaluated  for  the  special 
cases  considered,  A  summary  of  these  developments  follows  for  a 
special  form  of  the  F  function,  quadratic  in  y,  where 

.2 

F(y;y)  =  G(y)  +  y  Gi(y)  +  y 

Control -system  equation.  -  For  F  in  the  form  of  equa¬ 
tion  ( 44) , ^he  control- system  equation  (20)  becomes 


.2  ^  G(y) 

GgCy) 


(52) 


The  function  ^^(y)  does  not  affect  the  control  system  and  the 

control -system  equation  gives  an  explicit  expression  for  y  as  a 
function  of  y. 

Evaluation  of  integrals.  -  If  an  integral  of  a  func¬ 
tion  H(y,y)  is  any  one  of  the  criterions  to  he  considered,  it  can 
he  evaluated  as  follows: 


Thus  the  integrals  can  he  evaluated  without  solving  the  differen¬ 
tial  equation  (52) . 

Transient  behavior.  -  The  differential  equation  (52)  can  he 
easily  solved  as  follows: 


^0 


(54) 
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Degree  of  mninium.  -  Equation  (47)  was  derived  to  prove  the 
absolute  minimum  for  the  special  cases  considered  and  can  be  used 
to  evaluate  the  degree  of  the  minimum  found.  Thus, 


rit. 


H(y  ,  j*)  dt*  - 
0  UO 


J. 

H(y,y)  dt  =  GgCy)  Fi^Cy) 

*  -»r  L- 


fiti 


where 


H(y,y)  <it  is  to  be  a  minimum  or  maximum. 


Uo 


If  ^  ^  ^0  ^  ^  ^  ^  then  an  absolute  minimum  is 

obtained;  if  Cr2^y)  <  0  for  <  y  <  an  absolute  maximum  is 

obtained.  E<iuation  (55)  also  indicates  that  the  degree  of  minimum 
(or  maximiim)  varies  with  the  magnitude  of 

denominator  indicates  that  for  any  deviation  [y*-y|  in  y,  it  is 
better  to  err  on  the  side  of  a  larger  |;^  j. 


SUMMARY  OF  RESULTS 

A  rational  analytic  method  for  the  design  of  automatic  con¬ 
trol  systems  has  been  derived.  Criterions  on  the  behavior  of  the 
controlled  system  were  developed  in  the  form  of  certain  time  inte¬ 
grals.  When  any  of  these  arbitrary  but  physically  realizable  cri- 
"berions  were  used  as  a  starting  point,  those  equations  that  the 
control  system  must  follow  were  derived.  The  criterions  developed 
required  the  minimization  of  certain  time  integrals  using  the  cal¬ 
culus  of  variations .  The  method  gave  not  only  a  description  of 
the  behavior  of  the  controlled  system  but  also  gave  those  physicailly 
realizable  equations  that  the  control  system  can  follow  in  order 
to  realize  this  behavior.  General  results  were  shown  for  systems 
of  second  order  and  of  any  number  of  degrees  of  freedom. 


Lewis  Flight  Propulsion  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Cleveland,  Ohio,  October  11,  1950. 
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APPENDIX 

SECOND-ORDER  SYSTEMS  OF  SEVERAL  DEGREES  OF  FREEDOM 


It  is  beyond  the  scope  of  this  report  to  detail  the  cases  of 
higher-order  systems  and  those  having  several  independent  variables. 
The  general  equations  for  these  cases  will  be  developed  and  it  is 
expected  that  subsequent  reports  will  cover  their  applications  in 
detail . 


For  the  case  of  a  second-order  system  with  two  degrees  of 
freedom,  equation  (ll)  becomes 


nt 


F(y,y,y,z,z,*z’)  dt  =  minimum 


(Al) 


where  y  and  z  are  independent  f\mctions  of  time.  The  condition 
to  satisfy  equation  (Al)  is 


de 


n(t|]_+c5t 

.  II  It  ^  _  TI  It 

F(y+eSy,y+e&y,y-+€5y,z+£&z,z+€6z,z+£8z)  dt  =  0  (A2) 

Ut=0 


at 


£  =  0 


The  time  duration  of  the  integral  of  equation  (Al)  is  that  begin¬ 
ning  at  a  definite  time  (t=0)  but  not  ending  at  any  definite  time 

but  rather  along  some  curves:  y  =  f2_(t),  y  =  f2('t),  z  =  g]_(t), 
and  z  =  g2(t)  (fig*  H)  •  functions  By  and  Bz  are  arbi¬ 

trary  and  independent  functions  of  time. 

Performing  the  operation  indicated  by  equation  (A2)  gives 


(A3) 
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After  integration  by  parts. 


'^t3_ 

rv  2  1 

0 

By  dt  +  1  ^  ^  (^z^ 

UO  dt 

6z  dt  + 


(Fy  6y  +  Fy  ^ 


[f-  &z  +F-5Z  -  A(f.^.)6z] 


=  0 


As  before, 
condition  terms 
end  condition 


the  integrands  of  the  integrals  and  the  boundary- 
must  vanish  separately.  From  the  geometry  at  the 


(t=t2,  fig.  11) ^ 

=  1 

-  j(ti)] 

St^ 

&y(tj)  = 

[f2’(ti)  -  y(ti)] 

Btn 

V 

' 

Sz(ti)  = 

[gl'(ti)  -  z(ti)] 

Btj^ 

6z(t2)  = 

1  ^■*^1 

J 

The  three  conditions  from  equation  (A4)  then  become 
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F  _  _L  (p.)  +  (F-)  =  0 

at  ^  ^^2  ''  7* 


Fz  -  ^  (^z)  +  S  (^z)  =  0 

ax  ^54*^ 


F  -  ifj  +  i  ^  (Fy)  -  yFy-  -  iFJ  +  i  ^  (F”)  - 

fl'(tl)rFj.  -  ^  (Fy^  _  +  fE'(tl)  Fy(tl)  + 

L  J  t=t2^ 

&y(o)jF^  -  —  (F”)j^  ^  +  &y(o)  F”(o)  + 


>(A6) 


Si 


5z(0)  F*  -  ^  (F”)  +  5z(0)  F"(0)  =  0 

Jt=o 


The  first  two  of  equations  (AS)  are  the  differential  equa¬ 
tions  that  satisfy  the  original  criterions  of  equation  (Al).  The 
physical  solution  to  the  problem  is  the  pair  of  solutions  of  these 
equations  that  satisfy  the  boundary-condition  equation  of  equa¬ 
tion  (AS) .  If  the  first  of  equations  (AS)  is  multiplied  by  y 
and  the  second  by  z  and  the  equations  are  added,  an  exact  deriv¬ 
ative  is  formed  the  integral  of  ■which  is  as  follows; 

^  -  y^y  +  y  ^  <^y^  -  y^y  -  ^^Z  +  ^  ^  (FV)  -  =  c  (A7) 

A  study  of  the  boundary- condition  equation  of  equation  (AS)  then 
shows  that  the  physically  reasonable  boundary  conditions  should  be 
as  follows; 
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If 

F* 

^y 

If 

F’* 

'y 

If 

If 

F” 

Then 

(AS) 


problem  of  equation  (Al)  for  the  boundary  conditions  of  equa¬ 
tion  (A8)  is  as  follows: 

F  -  +  y  4  (Fv*)  -  ’yFy  -  ^  (^z)  -  *^Z  ^ 


dt  '  y 


dt 


and  either 


F 


y  dt  '  y 


(F,;)  +  ^  (F^o  =  0 


>  (A9) 


dt 


or 


The  meaning  of  the  houndary- condition  e<iuation  (A8)  is  to 
define  the  original  criterions  for  that  duration  during  vhich  the 
system  moves  from  one  essential  operating  level  to  another.  Thus^ 
if  all  conditions  of  equation  (A8)  must  hold,  the  system  goes  from 
one  definite  y,  y,  z,  and  z  to  any  other  definite  y,  y,  z, 
and  z.  The  first  differential  equation  of  equation  (A9)  vould  he 
of  third  order  and  the  second  or  third  equations  of  (A9)  vroiald  he 
of  fourth  order.  If  equations  (A9)  are  integrated  to  obtain  a 
pair  of  second-order  differential  equations  having  three  constants 
of  integration,  the  choice  of  these  three  constants  can  then  deter¬ 
mine  a  desired  end  point;  that  is,  the  values  of  y^  z,  and  z  at 
some  final  y. 

The  physical  solution  to  the  problem,  then,  is  the  pair  of 
second-order  differential  equations  that  are  solutions  of  equa¬ 
tion  (A9)  and  the  constants  of  integration  of  which  are  evaluated 
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so  that  the  system  goes  through  some  desired  end  point.  This  end 
point  must  he  such  as  to  allow  the  possibility  of  stability.  Such 
an  end  point  may  be  written  as  follows: 


when 


then 


y  =  ys 


y  =  0 

Z  =  Zg 

z  =  0 


^  (AlO) 

J 


vrhich  gives  three  conditions  for  the  evaluation  of  the  three  con¬ 
stants  of  integration.  A  stability  device  must  still  be  added  to 
the  system  so  that,  at  the  point  when  equation  (AlO)  holds. 


and 


*  *  A 

y  =  0 


z  =  0 


(All) 


Equations  (A9)  are  symmetric  in  the  variables  y  and  z, 
which  indicates  the  nature  of  the  extension  for  more  independent 
variables.  Thus,  if  a  third  independent  fimction  w  exists,  the 
original  criterion  would  be  written  as 


F(y,y/y,z,z/2f,w,w,w)  dt  =  minimum 


(A12) 


This  condition  is  satisfied  under  boundary  conditions  similar  to 
eq.uation  (A8),  where  two  additional  conditions  are  added  on  the 
variable  w  and  the  following  equations  describe  the  controlled 
system: 
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^  ^  "'y'  -  ■*'y  ’  H 


zF**  - 
z 


^v-  +  w  —  (F")  -  *^’  =  0 

w  dt  '  w  w 

and  any  two  of  the  following  three  equations; 

’'y  -  s  ^  '  “ 


^  "  dt 


F„  -  (F-)  +  ^  (F-)  =  0 

dt 


>(A13) 


The  three  equations  of  equation  (A13)  can  then  "be  integrated  to 
give  three  second-order  differential  equations  where  the  five  con¬ 
stants  of  integration  are  evaluated  so  that  the  system  goes  throu^ 
the  desired  value  of  y>  z,  z,  w^  w  for  some  final  value  of  y. 
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(a)  Speed  control.  (^)  Speed  overshoot. 


(c)  Temperature  overshoot  (d)  Compressor  surge, 

and  undershoot. 


(e)  Blow-out.  (f)  Compressor  line 

below  surge. 


-  Arbitrary  weighting  functions  for  various  control  criterions 
and  pertinent  engine  characteristics  Involved. 


Figure  2. 
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Temperature  integral 
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Temperature 


Figure  6.  -  Typical  transient  of  controlled  system  for  speed  control  with  temperature 

limiting  crlterlons  (case  B),  where  E  =  ^^2, 
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Figure  7.  -  Schematic  diagram  of  control  for  turbojet  engine  for 
minimum  rise  time  with  temperature-limiting  crlterlons  (case  C). 
Stability  device  must  be  added  to  this  system. 
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gure  8.  -  Various  control  parameters  as  function  of  A  for  minimum  rise  time  with 
temperature-limiting  crlterlons  (case  C),  when  accelerating  to  limiting  temperature. 
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